Single Molecule Photon Statistics from a Sequence of Laser Pulses 
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There are many ways of calculating photon statistics in quantum optics in general and single 
molecule spectroscopy in particular such as the generating function method, the quantum jump 
approach or time ordering methods. In this paper starting with the optical Bloch equation, within 
the paths interpretation of ZoUer, Marte and Walls we obtain the photon statistics from a sequence 
of laser pulses expressed by means of quantum trajectories. We find general expressions for Pn(t) 
- the probability of emitting n photons up to time t, discuss several consequences and show that 
the interpretation of the quantum trajectories (i) emphasizes contribution to the photon statistics 
of the coherence paths accumulated in the delay interval between the pulses and (ii) allows simple 
classification of the terms negligible under certain physical constraints . Applying this method to the 
concrete example of two square laser pulses we find the probabilities of emitting 0,1 and 2 photons, 
examine several limiting cases and investigate the upper and lower bounds of Po(i), and Piifi) 

for a sequence of two strong pulses in the limit of long measurement times. Implication to single 
molecule non-linear spectroscopy and theory of pairs of photons on demand are discussed briefly. 

PACS numbers: 42.50.-p, 42.50.Ar, 42.62.Fi, 82.53.Hn 



INTRODUCTION 

The interaction of matter with a sequence of laser 
pulses is a powerful tool frequently used for the investiga- 
tion of a wide variety of chemical, physical and biological 
systems [if. This field of research called non-linear 
spectroscopy uses clever design of laser pulses for the 
investigation of fast dynamics (e.g. pico - seconds) of 
ensembles of molecules in the condensed phase. Recently 
van Dijk et al Q reported the first experimental study 
of an ultra-fast pump-probe single molecule system. 
Unlike the previous approaches to such non-linear 
spectroscopy where only the ensemble average response 
to the external fields is resolved, the new approach yields 
direct information on single molecule dynamics, gained 
through the analysis of photon counting statistics 
In [sl we considered the non-linear spectroscopy for a 
single molecule undergoing stochastic spectral diffusion 
process. Here we neglect all dephasing and spectral 
diffusion effects, and concentrate on the effect of the 
laser field parameters on the photon statistics. 

Another related application is the generation of 
two indistinguishable photons using two short laser 
pulses interacting with a single molecule or atom 
Numerous applications for such photon sources have 
been proposed 0, 0, H, @| for the investigation of 
entangled states between identical photons and quantum 
properties of light, in the field of quantum informa- 
tion and quantum computation requiring consecutive 
photons to have identical wave packets. Usually in the 
mentioned experiments the single emitter is requested 
to supply one or two photons within as short as possible 
time interval. Although it is well established how to 
generate two photons from two ideal 7r-pulses if the delay 



interval between the pulses is very long, we can never 
produce two photons with probability equal 1, when the 
interaction time is finite. Thus, the information on the 
upper and lower bounds of the probabilities of emitting 
0, 1 and 2 photons as a function of the laser field pa- 
rameters, obtained in the manuscript, can be very useful. 

Although the th eory of single particle photon statistics 



ory 

is well established 10|, llll , it remained unapplied due to 



the absence of experimental ability to check the results. 
Recent experimental achievements 0, i, B m allowed 
the investigation of the interaction of a single quantum 
system with an external laser field inspiring further devel- 
opment of theoretical methods [3, [l3, 3, 



HE, EE 



Today there are several approaches to photon counting 
statistics such as generating function method [l^ or 
quantum jump approach [20| suitable for analytical 
predictions and numerical calculations. In this paper 
we follow the path interpretation approach of MoUow 
and ZoUer, Marte and Walls lllj of the optical 



Bloch equations [2lj, and show that this method is 



very useful for the analysis of single molecule non-linear 
spectroscopy. 

In what follows we consider a two level molecule 
interacting with two laser pulses and obtain general 
expressions for Pn{t) - the probability of emitting n 
photons in interval (0, t) by means of quantum trajec- 
tories. We discuss the influence of the coherence on the 
photon statistics. Also the explicit calculation of Pq, Pi 
and P2 - the probabilities of emitting 0,1 and 2 photons 
in the limit of long measurement times {t — > 00) is 
investigated in detail using the example of two identical 
square pulses. Some technical details skipped in the text 
are given in Appendixes A, B and C. 
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where the Rabi frequency is Q 



id 



Eg and d„e is 



OPTICAL BLOCK EQUATIONS 

Interaction of an atom or a molecule with a radiation 
field is described by the optical Bloch equations under 
well established conditions [2l[, and we remind the 
reader some of the basic assumptions. First (i) the laser 
field is intense, so that it can be modeled classically. 
Here the external electric field is E(t) = Eo/(t), where 
the amplitude Eq is independent of time. (ii) The 
electronic states of the single emitter are modeled based 
on the two level system approximation. This assumption 
is excellent when the laser is resonating with a particular 
absorption frequency of the molecule, the latter being 
well separated from other natural frequencies of the 
emitter. Most single molecules have a triplet state, 
however the life time of the triplet is much longer than 
the time scales under consideration in this manuscript, 
and it can be neglected, (iii) The spontaneous emission 
process is described by the Markovian approximation, 
where the inverse life time of the excited state is F. (iv) 
We neglect thermal dephasing, spectral diffusion and 
interaction of the emitter with a thermal bath, which 
was partially treated in [3]. (v) Finally, we will assume 
that the dipole moment of the excited and ground state 
of the single emitter is zero, so that only the transition 
dipole moment of the particle is important. Assumptions 
(i,ii,iii,iv) are physical assumptions which are justified 
in many single molecule experiments at least at low 
temperatures [l^ [2^ . and condition (v) is not limiting 
since our technique could be modified in principle to the 
case where excited and ground states of the molecule 
have a dipole. 

The two level system is described by a vector composed 
of the density matrix elements: a = (ctog, fgg, tige, Ccg)"^- 
Here CToo and dgg represent the populations of the excited 
and ground states respectively and CgcCTog describe the 
coherences, namely the off diagonal matrix elements of 
the density matrix, and obey cr*g = age- The optical 
Bloch equation is [21] 
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the transition dipole moment of the two level system. 
The operator F Eq. ([3]) describes direct transition from 
the excited to the ground state, and hence is associated 
with the spontaneous emission of a single photon. 

The optical Bloch equation Eq. ^ does not yield a di- 
rect method for calculating the probabilities of the num- 
ber of emitted photons. However starting with 
an interpretation of the optical Bloch formalism yields a 
tool for the calculation of photon statistics, based on the 
n-photon- propagators (see details below) . The formal so- 
lution to Eq. ( HI) may be given by the infinite iterative 
expansion in F 



(7(4) = e(t, O)a(o) + / dti5(i,<i)Fg(ti,0)a(o) + 



d<2 / dti^(t,<2)fSfe,ii)fe(ti,0)a(o)- 
Jo 



(4) 



where f7(o) is the initial condition, and the Green function 
describing the evolution of the system in the absence of 
spontaneous transitions into the ground state (i.e. with- 
out F ) is 



e(t, i') = rexp 



L(ti)dti 



(5) 



where T is the time ordering operator. The first term 
in the expansion Eq. ([4]) does not include F at all, and 
hence describes the process where no photons are emit- 
ted, the second term includes F just once and describes 
the process where one photon is emitted etc. It is there- 
fore useful to define the conditional state a'^^ , where n is 
an index for the number of photons emitted in the time 
interval (0, i). Then by definition 



(") 

(t,0)'^(0)' 



(6) 



where the n-photon-propagator [l5| is 



(n) 
(*,«') 



Air. 



dti g(i,t„)f •••fe(ii,i')- (7) 



The physical origin of the n-photon-propagator defined 
by Eq. ([7]) is simple and intuitive: the system evolves 
interacting with the laser field without photon emissions 
until time t\, it then emits a single photon and contin- 
ues the evolution without emissions until time ti when 
it emits the second photons and so on. At this point it 
is convenient to choose a four-dimensional orthonormal 
basis to work with: |e) = (1,0,0,0)'^, \g) = (0,1,0,0)^, 
|c) = (0,0,1,0)^ and |c*) = (0,0,0,1)'^. According to 
the matrix form of the Bloch equation Eq. ([T]) the first 
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two vectors correspond to pure excited and ground states 
respectively. The last two vectors, however, do not rep- 
resent any real physical state and should be simply con- 
sidered as convenient mathematical way to include all 
possible quantum paths going through superposition of 
the pure physical states |e) and |(?), thus representing the 
contribution of the coherence effect. Using this notation 
the main equation for calculating the probability of n 
emission events up to time t is [23| 



Pnit) - ((e| + (g|)cT[,"))) = ((e| + {gDUH'^^W^o)) ■ (8) 
For example the probability of emitting zero photons is 



Po{t)= ^(z|e(i,0)|a(o)), 



(9) 



Eq. ((T4)) means that the propagator corresponding to 
n emission events in {t',t) can be decomposed into a 
emission events in (f , to) and n — a emission events in 
{ta,t). The extension to more than one time point such 
as ta is trivial and leads to summation over all possible 
permutations of the n photons propagators resulting in 
n emission events. 

Turning back to the Eq. ([H]) for P„(i) and inserting the 
closure relation 



(15) 



we find 



and the probability of emitting a single photon is 

Piit)=y2mf dtig(i,ti)f^(ii, 0)1(7(0)). (10) 

._ JO 



Consider a laser field interacting with the molecule in 
the time interval (t' , t) and choose a fixed point ta inside 
this interval. Such a partitioning of the time axis is use- 
ful for the analysis of sequence of pulses investigated in 
the following section, when we distinguish between time 
intervals where the laser is turned on and off. First, let's 
split the integration over t„ in Eq. ([7]) into two parts: 



dtn 



dtn ■ 



(11) 



Using the fact that in the first interval {t',ta) tn < ta 
and replacing the Green function Q{t, t„) by the product 
Q{t,ta)Q{ta,tn) One easily finds 



dtr, 



At^g{t, tn)t ■ ■ ■ vg{t^X) = 

(12) 

Now, left with the integral over the second range {ta,t), 
we repeat exactly the same procedure as we did with the 
initial expression, but this time we split the integration 
over tn-i into (<', ta) and {ta, tn). Similarly, using t„_i < 
ta and replacing ^?(i„, i„_i) by C?(i„, ia)^(ia, ^n-i) inside 
the first interval we find 



dtn r Atn-i--- f'dtig{t,tn)t---tg{ti,t')^ 

ta JV Jt' 



(13) 



Repeating this algorithm n times it is easy to prove that 

n 

j{n) _ jj{a) jj{n-a) ^^^^ 



(*,*') ~ {t,ta)'-^{ta,t') ■ 



a=0 



i—e,g j—e,g,c.c^ a— 

(16) 

Eq. (|16p describes the summation over all possible paths 
resulting in n emission events and suggests the following 
classification : the paths going through the pure states 
\j) — \g) may be identified as semiclassical, whereas 
the paths going through the states \j) ~ |c), |c*) describe 
the contribution of the coherence. 



TWO PULSES 

Now we focus on the case of two laser pulses separated 
by a window A in which the laser is turned off. The initial 
time is t = 0, the time ti is the moment when the first 
pulse is switched off. The amplitude of the external field 
remains equal zero f{t) = for the delay period ti < t < 
ti + A. At time ^2 = + A the laser is turned on again, 
and then again turned off for t > t^ {f{t) = for t > 
ta). Schematically the sequence is represented in Fig. 1 
for square pulses, however we emphasize that the results 
obtained in this section are valid for pulses of any shape. 
Our goal is the derivation of general expressions for P„ (t) 
from two pulses in the limit of the long measurement time 
t oo when we know that eventually the system is in 
the ground state. We assume that the molecule is always 
in the ground state at the beginning of the experiment. 
If we divide the time axis into four distinct intervals : 
two intervals when the laser is turned on and two others 
when the laser is turned off, the most general expression 
for ^^("0-) following from the extension of Eq. is 

ttM _ 7-r(n-"-/3-7)rr(7) TiiP) Tj{a) /i 

where the superscripts a, /3 and 7 are all non negative 
integer values leading to n photons (i.e. n — a — /3 — 7 > 
0). The Einstein's summation rule from to n must 
be applied to every superscript appearing twice. Inside 
time intervals {t, t^) and (^2, ^i), when the laser is turned 
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ofF, the Rabi frequency is equal zero = 0, and the 
calculation of the Green function Q{t, t') Eq. (O becomes 
nearly trivial. For the delay interval A we find only two 
non-zero n-photon-propagators: 



Tj{0) 



-FA 










g(i:i^o-r/2)A 








-(iwo+r/2)A 



U, 



(1) 

(ti+A,ii 



/ 0\ 

l-e-^^ 





(18) 



(19) 



and U, 



(n) 

(ti+A,ti) 



for n > 1. This result is definitely 



expected , since if nothing excites the molecule, there is 
no chance to get more than a single photon. The matrix 
representation of the propagators in the interval {t > 
ts) when i — > cx), are found by taking the limit A — s- 
oo of Eqs. (fT8|l .(|19 p . Now inserting the closure relation 
Eq. (fT5]l between each two propagators of Eq. (fT7|) and 
using the just obtained matrix elements Eqs. P^ . p^ 
we find : 



where 



Q=0 



lim Pn{t) = P^^^ + e-^^(e*^"«A^°'^ + C.C.), (20) 



{(5l<r;i5>(5l<o)l5) + e-'^^(<7l<r;|e)(e|<„)|g>} + ^(1 " e-^^)(ffl<rr'^l5)(e|<lo)li/) 



n-l n-2 

+ E {(^l<r)"^l^')(5l<o)l5> + e-'^^(e|<-:)-^)|e)(e|<„^|5)} + E(l - e-^^)(e|C/£-^)-^^|5)(e|C/£;o)|5> 



and 



- E(5i<:.Si^)(^i<o)i5> + E(^i<rr'^i^)(^i<o)i5>. 



(21) 
(22) 



0=0 



a=0 



Eqs. (201), (ED) and summarize all possible paths 
resulting in n photon emission events and allow sim- 
ple identification of negligible terms, when particular 
physical constraints are taken into account. It is easy 
to see, that the first two terms of P^''^ Eq. ([2T|) (those 
with J2a=o) describe processes where all n photons 
are emitted during the pulses and none in the delay 
interval or after the second pulse. The third term 
of this expression represents the processes, where a 
single photon is emitted in the delay period (with 
probability [l — e"'^'"]) and n-l photons during the 
pulses events. Similarly, the next two terms originate 
from the processes, where a single photon is emitted 
after the second pulse and zero photons in the delay 
interval, and finally, the last term describes situations, 
where one photon is emitted in the delay interval and 
another after the second pulse. This interpretation may 
be used to simplify the calculations, as for instance in 
the case of the short pulses considered below, where we 
neglect the trajectories with photons emitted during the 



pulse events. 

Ahhough Eqs. ([201), (ED) and ^ are very general, 
they already contain interesting physical information. 
First of all, we pay attention to the fact, that the 
coherence terms A^°^, describing the processes where 
the molecule is left in the superposition of the pure 
states at the end of the first pulse (i.e. the paths going 
trough the states |c) and |c*)), never include trajectories 
where a photon is emitted within the delay interval A. 
Mathematically this follows from Eqs. ([T8|) and (fT9|l . and 
physically it makes sense, because spontaneous collapse 
into the ground state destroys the coherence. Secondly, 
we see, that the coherence terms are multiplied by the 
exponentially decaying factor e"^'^^^, responsible for 
the dephasing effect, and oscillate in A with orbital 
frequency ujq (see the e*^'^° term in Eq. (|20|) ). In 
optics luq is much larger than the inverse of r - the 
minimum time resolution of the measurement device 
: tluq 3> 1. Therefore, in order to match our results 
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for the probability of emitting n photons to those 
observed by an experimentahst, it is essential to treat 
the coherence terms as stochastic variables - i.e. it is 
reasonable to replace them with their time average, 
which is equal zero. However, it should not be forgotten 
that : (i) in the limit A ^ 0, when the pulses are 
attached together, the coherence contribution A'^°^ 
becomes non-oscillating and non- negligible part of P„, 
and (ii) in non-optical microwave experiments, where 
the absorption frequency is comparable with the time 
resolution of the measurement device the influence 
of the coherence trajectories is important. 

It is possible to derive another useful expression for the 
probability of emitting n photons from two pulses. First 



we note, that according to Eqs. ([8]), (jlSp and (fTQ]) the 
probability of emitting n photons from any single pulse 
or sequence of pulses of total length T in the limit of 
infinitely long measurement time t — s- cx) is 

^» = (5|f/(S)l5> + (e|f/fT,of^l<7>, (23) 

(where for n=0 the second term (e|C/^^j, I5) = 0). From 
the physical point of view the second term of Eq. ((23|) 
expresses the fact, that the molecule left in the pure 
excited state eventually decays to the ground state by 
spontaneous photon emission. Simple rearrangement of 
Eqs. (|2ip and ([22|l . with details given in Appendix A, 
results in : 



Eph ph 
n — OL a 

a=0 



-Ar 



Q=0 



„A(r/2+ic^o) 



Coh 



c.c. 



(24) 



\g) is the probability 



HerePi^ = (5|C/£!o)|ff) + (e|C/(';;o) 
of emitting n photons only from the first pulse [24 
similarly P^^^ = {9\u'i^lt^)\9) + {e.\u\^~^^^\g) designates the 
probability of emitting n photons only from the second 
pulse, and 



0=0 



(25) 

is the probability of emitting n photons from the two 
pulses produced one immediately after another (i.e. with 
zero delay). This formulation of P„ Eq. (p4|) shows, that 
the first term J22=o Pn-aPa represents the sum of all 
possible ways of emitting n photons from the both pulses, 
as if the consequences of the interaction of the molecule 
with the first pulse had no influence on the state of the 
system at the beginning of the second pulse, i.e. like if 
the treatment of each pulse could be done independently. 
Nevertheless, since such an influence exists, it is reason- 
able to define the rest of the terms on the righthand side 
of Eq. p4)) as a correlation: 



C{/^)^Pn-YPto.Pt 



-,-Ar 



-I- 



(e 



A(r/2-|-iwo) 



(26) 



Note that Eq. (|24p makes perfect physical sense in the 
limits A — > cx) and A — > where we find trivially ex- 
pected results. In the first case only the first term on the 



righthand side of Eq. (p4)l survives - i.e. this limit de- 
scribes the situation where the interaction of the molecule 
with the first pulse indeed has no influence on the interac- 
tion of the molecule with the second pulse, since all coher- 
ence effects have enough time to decay completely. And 
the second limit gives P„ — Pn^^^ ■ We emphasize, that 
in the flrst case, once the probabilities of emitting n pho- 
tons from each single pulse are known, the efforts needed 
for the calculations are considerably reduced. However, 
care must be taken while using Eq. ([M)) for the calcula- 
tion of the second limit A — > 0, since the continuity of 
the laser's phase plays important role, as demonstrated 
on the example of two square pulses in the subsequent 
section. 



EXAMPLE : TWO SQUARE PULSES 

In this section we apply our general results to the con- 
crete example of two identical square laser pulses. Con- 
sider the sequence : 



/ J n—a a ' 



fit) 





cos [uJhit - t2) 





< i < ti 

ti <t <t2 
t2 <t <t3 

h < t 



(27) 



where ti = t^—t2 = T - is the pulse's duration, t2—ti = A 
- is the delay period between the pulses and (j)i and 02 
are the initial phases of each pulse. For simplicity we as- 
sume, that the laser frequency lol — namely we con- 
sider the case of zero detuning, and also the initial phase 
of the first pulse is zero <j)i = Q. The time dependence 
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t 























T 



to t 



(a) 



(b) 



(c) 



(d) 



t (Time) 



and A'^"^ Eqs. dH]) , ([221) by replacing all n-photon- 
propagators Eq. ([7]) by - the n-photon-propagators 

calculated within RWA (see the derivation of Eq. ((7T|) 
in Appendix B) . All mathematical manipulations and 
calculations were obtained with the help of Mathematica 
5.0. 

First we consider in detail the probability of emitting 
zero photons. For the calculation of Po{t) we need the 
zero photon propagator. This is the simplest case, since 
there is only one possible permutation. According to 
Eq. ini) we have : 



(0) 
t,0) 



(0) 



(0) 



(*,t3) (t3,*2) (*2,tl) (tl,0) 



(29) 



FIG. 1: Two square laser pulses as modeled in Eq. (|27[) . T is 
the length of a single pulse and A is the delay interval between 
the pulses. The dashed arrows show four cases out of eight 
possibilities of emission of two photons. For example: (a) two 
photons are emitted within the duration of the first pulse. 



Inserting the closure relation Eq. (|T5|) and applying RWA 
lead to : 



„CIa 



r(0) 



(t3,t.)ie)(e|'^(V;.o)ig^ 



of the Rabi frequency fi, describing the interaction, is 
schematically illustrated in Fig. 1. The propagators act- 
ing on the molecule during the square pulses are found 
using the Rotating Wave Approximation (RWA). In Ap- 
pendixes B and C we show, that within RWA Eq. (|20|) 
can be rewritten in the form 



+(gi<.,)ig)(gic/(:,:o)ig)' 



(0) 



«o^°' = (g|<..)|c)(c|C/(^,7,o)lg) 



~r(0) 



(30) 



(31) 



where the definition of p„ and a„° follows from 



Calculating the matrix elements of Eqs. (l30|) and ([3T|) 
we find the following explicit expression for Pq ~ 
limt^oo Po{t)- 



i^-T-A 



(1 - 4172 



(1-4172)' 



1 - 2Vl^) cosh 



/ VI - 4172 \ / ^1 _ 4f72 \ 

I T J + Vl^4f72 sinh ^-^ T J 



8172 



fl-4172) 



.e--/^smh^h^lHi^r 



yi^l^coshKlHi^T 



/ VI - 4172 N 
sinh ■ T 



COS[CJO (A + r)-(^2], 
(32) 



where we set F = 1 for simplicity. The last term, exhibiting oscillations due to the cos[wo (^ + ?") ~ 02]i 
results from the quantum paths going through the |c) and |c*), thus representing the coherence effect. Of 
course, when T = or 17 = 0, Pq = 1 since no photons are emitted, and if T — > cxd, Pq = Q since many photons 
are emitted. Similar calculations were made also for Pi and P2 - see Eqs. (|8T|) . (j88|) for the final results in Appendix C. 



For very intense laser fields, when the Rabi frequency is much larger than the inverse life time of the excited state. 
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taking the limit Q ^ 1 of Eq. (|32p we obtain : 



And using Eqs. ((8T|l . 



lim Po ~ e-^ |e-^ sin^ (^^^ + cos^ (^^^ - ig^^/^ sin^ (r!T) cos[^o (A + T) - (/.^ij 



(33) 



lim Pi = <j 4sin2(rjr) + 2r 

03>1 8 



1 - sin^ f ^) (1 - e-^) 



Tsin^inr) (l + e 



-(T+A/2) 

+ (T + 2) sin^nT) cos[cjo (A + T) - ^2] 



(34) 



lim P2 — e 



sm 



4 t2 



V 2 y 64 



cos(fir) 



[l~e-^) + j[cos''{nT) + 2] 



g-T -(T+A/2) 

(cos2(nT) + 4) (1 + e""^) r(T + 4) siii^{nT) cos[wo (A + T) 

32 16 



02]. (35) 



Finally, we would like to investigate the limiting be- 
havior of Po, Pi and P2 within the strong fields approx- 
imation in the case of long A ^ 00 and short A — + 



delay intervals. As shown in [2J,|25|, the probabilities of 
emission 0, 1 and 2 photons from a single square pulse of 
length T (see Eq. ((23l) ) are given by 



lim PJ. 

S2>1 



lim Pr 







lim P'^ 



-T/2 2 

e ' cos 



lim P/ 



-T/2 



and 



lim Pi 



-T/2 



[4 2r - (4 -I- T) cos(17T)] 



(36) 
(37) 



64 



-T [4r 16 (8 + T) cos{nT)] . (38) 



Taking the limit A > 1 of Eq. 



we find 



I\2 



lim Po = (PoO 
n>i,A>i 



-T 4 

e cos 



2 



^ , (39) 



which is equal precisely to the product of the probabilities 
of emitting zero photons from two single square pulses 
Eq. ([36|) and completely agrees with Eq. ([24]) . Now using 
Eqs. dsn, ([321), dSSl) and the fact, that the two pulses are 
identical, we can easily obtain the limit A — )■ 00 of Pi 
and Pj: 

lim Pi = 2P/Pi = 

n>i,A»i 



3T)cos^ — -2(4 + T)cosM — 



^7^ 



(40) 



and 



.-T 



n>M:A»i = ^2 + (AO = [24 + r(40 + 9T)- 



(-32 T^) cos(rJT) (8 T(8 T)) cos(2rjr)] . (41) 



Considering the opposite limit A — > we remind, that 
the contribution of the coherence paths going through the 
states |c) and |c*) must not be neglected. Moreover, it is 
essential to take into account, that two attached square 
pulses of the same Rabi frequency are equal to a single 
long square pulse, only if 4)2 = (t>i + wo(A + T) - i.e. the 
pulse is continuous. Assuming for simplicity, that this is 
the case, when A ^ from Eq. ((32|) we find for Pq : 



hm Po 

n»i,A<i 



(42) 



which once again agrees with Eq. ([24)) . since it is exactly 
the result of replacing T with 2T in Eq. (|36p. Similarly 
the limits A ^ of Pi and P2 may be obtained by 
replacing T with 2T in Eqs. ^ and ([SS]) . 



In Fig. 2, neglecting the fast oscillating coherence 
paths, we plotted the semiclassical terms Pq^'^, Pi^'^ and 
Pj'** for the relatively long A = 3 and short A = 0.5 
delay intervals for the case of strong laser field fl = 10. 
Comparing the graphs one may see, that the dependence 
of pf^'^ on A is visibly weaker than those of p^^'^ and P2^^ 
(we explain this effect later - see the discussion below Ta- 
ble 1). When FT > 1, we expect that : (i) Po''',_pf''' and 
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FT 




0.6 



FT 



0.8 



0.6 




FT 



n 




FIG. 2; The semiclassical parts of probabilities of zero, one and two (from left to right) photon emissions from the two laser 
pulses as a function of T - the duration of a pulse with fl — lOF. The solid line is the exact result Eqs. (|32|l . (ISlfl and (I88p 
and the dashed curve show the approximation in the limit of strong fields Eqs. p3|) . (|34|l and (|35|l . The upper row (a, b, c) 
illustrates the result for A = 3 and the lower row (e, f, d) illustrates the same probabilities for A — 0.5. 



P2^^ are all small, since many photons are expected to 
be emitted during the pulses, and (ii) independent of A, 
since the contribution of the photons emitted during the 
pulse event is much larger than the contribution of the 
photons emitted in the delay interval. Such a behavior 



is clearly seen for P2 (compare Figs. 2.c and 2.f ) where 
the difference between the case A = 3 and A = 0.5 is 
stronger for short T. Below we prove this in a Poissonian 
limit. 



Strong and Short Pulses 



following matrix representation : 



The sequence of two very short and strong pulses is 
important due to its numerous practical applications. 
Mathematically we define this limit as T ^ 0, ^ oo 
in such a way, that the product fiT stays of the order 
of unity flT ~ 1 (automatically leading also to 3> F), 
otherwise the molecule will never reach the excited state, 
and the probability to obtain non-zero results for Pi and 
P2 will be negligible. As a consequence, the spontaneous 
emission process during the pulse event may be neglected, 
and therefore, it is reasonable to approximate the behav- 
ior of the system by simple Shrodinger evolution with 
well-known Rabi oscillations. In this limit the photons 
can be emitted only in the delay interval or after the 
second pulse, while the only non-zero propagator acting 
on the molecule during the pulses within RWA has the 



lim U^°'^ 

O — ^CJO 



/ C0S2 ^ 

2 nT 



sm „ 

■ sin OT 



V 



sin2 ^n: 

cos^^f 

„■ sin OT 

' 2 
• sin QT 



^ sin OT jsinOT 
■ sin f?T ■ sin f2T 



COS 

sin 



2 

2 SIT 



sm 



^# cos^t; 

(43) 

which is independent of F. Note, that since this zero- 
photon-propagator Eq. describes the conservative 
evolution of the system, the transformation is unitary, 
all the elements exhibit Rabi oscillations, and symmetry 
and reversibility of the matrix elements are found: 



(e\U, 



(0) 

(T,0)l 



(T,0)ly/' 



{g\u[%^\e) etc. 



Now we consider Pq, Pi and P2 in the limit of the 
short and strong pulses and demonstrate how this phys- 
ical constrain can help in reducing the number of paths 
appearing in Eqs. (P0)) . (PT|) . (|22p . The exact expressions 
for C/i^n\ and J7,^f L , according to Eq. pT|) consist of the 



(t,o) 



(t,0)' 
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sum of 4 and 10 terms respectively - see Eqs. ((78| . (1861) in 
Appendix C. After neglecting the trajectories where pho- 
tons are emitted within the pulse events (since T ^ 0)) 
we have : 



tAO) tAO) ttW rr(0) 
■ ^ (t.tsr it3,t2r it2.ti)^ (tuO) 



(44) 



and 



This is one example where the formulation of photon 
statistics based on quantum trajectories is very con- 
venient, since we can identify the underlying physical 
processes and make approximations. Inserting the clo- 
sure relation Eq. between every two propagators of 
Eqs. (UH), ([Is]) and using the matrix elements of the zero- 
photon-propagator Eq. (j43p . we obtain the leading semi- 
classical and coherence terms of Pq , Pi and P2 in the limit 
of the short and strong pulses. The results are summa- 
rized in Table 1 below and they are valid for any sequence 
of short pulses. 



lim 

■*o,n- 



Am 



(0) 



(1) 



(0) 



(tM) (HM) (t2M) (tifl) 



(45) 



n 


„Cla 







(glC) Is) (gl<,0) Is) + ^-^''^?>KL) l«) (^l<,o) Is) 


(glC)l^)(^l<,o)lg) 


1 


(glC)|g)(^l^(".Ulg> + (e|?^(^"j,.)lg)(gl^(*'!o)lg) 


(^l^feU)l^)(^I^S!o)lg) 


2 


(^l<..)|g)(^l<o)lg) (1-^-^^) 






Table 1: Photon statistics for short pulses. 

Note, that the coherence terms of P2 vanish, since emission of a photon in the delay interval, necessary for emitting 
two photons in the limit of very short pulses, destroys the coherence . Using the symmetry and reversibility of the 
zero-photon-propagator matrix elements Eq. it is easy to show that a^"^ = —a2°^ and pp'^ +pf + P2^^ = 1, 

i.e. the semiclassical paths conserve probability. Finally, we bring attention to the fact, that the semiclassical paths 
of Pi do not depend neither on the spontaneous emission rate F not on the delay interval duration A (see Fig. 2.b 
and Fig. 2.e). Comparing the two non-negligible trajectories of p^'' with Eq. ([H]), we see that they correspond to 
the first term of the righthand side of this equation - i.e to the product of probabilities of emitting and 1 photons, 
related to each one of the pulses independently. 

Using the matrix elements Eq. (|43)) and Table 1, after some algebra we obtain explicitly : 

T^^fUoo ^0 = (¥) + (?) - r (T + A) - ^2] (46) 

which is the ^ 00, T ^ hmit of Eq. ((5^ . 

hm Pi = \ sin^ {VlT) [l + e-^'^ cos [wq (T 4- A) - ^a] | (47) 



and 



l-e-^)sinM:^) (48) 



hm Po 



It is easy to see, that when T ^ the Eqs. ([34]), ([35 
reduce to Eqs. ([Tfl) . (H5)) as expected. 



For a TT-pulsc defined by VlT = tt + 27rn 2J, |25| , where 
n is a non negative integer, in the strong field limit we 
obtain : 



lim Pq ~ e 

T^O,n-*oo 



-A 



(49) 



This behavior may be easily understood for : substituting 
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riT = n into Eq. (gS]) we have 



P2 Eqs. mM 



we find 



lim 

nT=7r 



u: 



(0) 

(T,0) 



/ 1 \ 

10 

1 

V 1 / 



(50) 



lim Po = e-^^ 



lim Pi = 2r2^re^ 



(53) 



(54) 



The physical meaning of this propagator follows directly 
from its matrix representation : as well-known, the ideal 
TT-pulse simply switches the state of the molecule from 
the excited to the ground state and vice versa. Thus, 
the first TT-pulse of the sequence pumps the molecule 
from the ground state to the excited state. If the de- 
lay between the pulses A is long, the molecule will 
emit a photon before the arrival of the second pulse, 
and then limT^o.n^oo Po = 0- Contrary, if A <C 1, 
the second 7r-pulse pushes the molecule back to the 
ground state before the emission of a photon, and then 
limT^o.a^oo Po = 1- We also notice, that the interfer- 
ence cos [luq (P + A) — term vanishes for the 7r-pulses, 
since according to Eq. ((50|) the off-diagonal matrix ele- 
ments giving raise to this term are equal zero. On the 
opposite, for a 7r/2-pulse, defined by flT = 7r/2 + 2nn, 
the influence of the coherence on the photon statistics 
generally does not vanish: 

lim Po - T (e"'^ + 1 - Se-'^/^ cos [wo(A -I- P)]| (51) 



Once again we remind, that in optics in many cases 
the ideal tt and 7r/2-pulses are considered where the 
interaction time FT —> 0, and then e~'^ = 1, but ujqT 
is not a small number, especially because we work 
under the assumption Q <ti loq, which is essential for 
the two level model approximation of the molecule 
and for the assumption, that the spontaneous emission 
rate F is not effected by the presence of the laser field 21 1 . 



fim P2 

n^O.T^oo 



(2i72r)2 



lim 

n^o,T- 



P„ 



(2f72r)" 



(55) 



(56) 



Indeed one can show, that Eq. (|52p and the limiting 
behavior Eqs. (|53ll56| of the exact results are identical. 
Note, that Eq. ([56]) corresponds to Poissonian statistics. 
This behavior originates from the fact, that because of 
the long pulses duration and weak laser field, the leading 
terms of P„ are those, where photon emissions are well 
separated one from another on the time axis, and there- 
fore photon statistics is described by nearly uncorrelated 
emission events. 



The upper and lower bounds for strong pulses 

Finally we investigate the upper and lower bounds of 
Po, Pi and P2 within the strong field approximation. Cal- 
culating the first and the second order partial derivatives 
with respect to T, we find the extremum of Eqs. ((33|) , ((34|) 
and (|35|) . and neglecting the ultra- fast oscillating coher- 
ence terms, obtain the bounds of semiclassical parts of 
Po, Pi and P2 : 



-(T+A) < 



„Cla 

Po 



< e 



(57) 



Weak and Long Pulses 

Here we consider the case of very long and weak pulses. 
In this limit the delay interval has a negligible effect on 
P„. In this limit we expect, that according to Eq. ([^1]) 
the probability of emitting n photons from two separated 
pulses may be approximated by P^^^^ - the probability 
of emitting n photons from two attached pulses with zero 
delay : 

n n — 1 

^linr^P„ ^ Pi^'^ = + ^(^1^(27)15) 

Q=0 a=0 

(52) 

Taking the limit — > 0, T — > cx) in such a way that 
il^T remains finite of the exact solution for Pq, Pi and 



^re-^<p^'^<e-^(l-e-^ + ^r). (59) 

where once again F = 1. 

The origin of Eqs. (HT]) , dSS]) and dM]), although non- 
trivial for a finite T, can be easily understood in the limit 
of the short pulses T ^ 0. The upper bound of Po is 
obvious, since if the interaction time is zero, no photons 
will be emitted for sure. Further, since we neglect the 
probability of emitting photons during the pulses, Po'*^ 
may only be decreased by the probability of not emitting 
a photon during the delay interval - e~^. Considering 
p^i^ we see, that Eq. (|T7)) reaches the maximum value 
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- 5) when QT = 7r/2. Hence, the maximization of pf^^ 
corresponding to the interaction with a sequence of two 
pulses is achieved by applying two ideal 7r/2-pulses. 
Finally, from Eq. ((48|) for limr^o.n-^oo P2 follows, that 
maximum of this expression is found for QT ~ tt, and 
hence the optimization of P2 is achieved by two ideal 
TT-pulses. 

Learning from this simple example, although not 
rigorously, we extend it to the conclusion, that the 
maximum of (n > 1) for any fixed interaction time 
is optimized by producing n equally separated 7r-pulses. 
This statement follows from the following argument : as 
far as we work under assumption, that the laser field 
does not effect the spontaneous emission process, the 
maximization of P„ in any limited time interval may be 
achieved by minimization of induced emission, which is 
guaranteed by the strong and short ideal 7r-pulses better 
than any by any others. 

In Fig. 3 we show the maximum of from a se- 
quence of two equal square pulses as a function of the 
interaction's strength and pulse's duration T for three 
fixed values of the total interaction time 2T + A . The 
curves were obtained using the extremum conditions of 
the exact expression for P2 Eq. l|88|) . From Fig. 3 we 
see, that as the interaction time FT becomes shorter, the 
delay period A longer and the Rabi frequency Q larger, 
the probability of emitting 2 photons is getting close to 1. 
Although it becomes equal exactly 1 only for two ideal tt- 
pulses separated by infinite delay, the graph shows, that 
starting from some range of parameters the increasing of 
P2 slows down, so that further increasing does not 
contribute much. Finally, we note that for short delay 
intervals A < 1 the maximum of is much less than 
1 as expected. 



SUMMARY 

We obtained general expressions for the probability of 
n photon emission events for a two level system interact- 
ing with two laser pulses separated by a delay interval A. 
The photon statistic was represented as summation over 
quantum trajectories, which allowed simple intuitive 
physical interpretation of the final results Eqs. ([20|l . ([2T|) . 
([22| . In particular, the contribution of the coherence 
effect, resulting from the quantum trajectories going 
through the superposition of the pure states at the end 
of the first pulse, was discussed. Although in optics it 
might be difficult to detect this effect experimentally, 
since the coherence paths oscillate in A with extremely 
large molecule absorption frequency luq, nevertheless 
in microwave spectroscopy, dealing with lower range of 
absorption frequencies, the coherence effect is important 




FIG. 3: The maximum of the probability of emission of two 
photons in a two-pulse laser field. The ts of the end time of 
second pulse is fixed at 0.5 (solid curve), 2.0 (dashed curve), 
4.5 (dot-dashed curve), respectively. The delay time between 
two pulses, A — — 2T. The star gives the asymptotic 
behavior of Eq. (|88[) in the limit of strong fields (57 =50) and 
short pulse with fiT = tt, F = 1. 



8]. In addition, the correlation function C(A) Eq. 
was suggested as a measure of the photon statistics 
deviation from a treatment where the sequence of 
pulses is considered as if the pulses were independent. 
This correlation might be a useful tool to quantify the 
coherence and "memory" of single molecules, atoms 
or quantum dots through the measured photon statistics. 

The application of our general results was demon- 
strated on detailed calculation of Pq, Pi and P2 - the 
probabilities of emitting 0, 1 and 2 photons from the 
sequence of two square laser pulses Eqs. (|32p . ([8T|) and 
(|88p in the limit of long measurement times t 00. The 
physical interpretation of the quantum paths was shown 
to be useful in reducing the complexity of calculations 
in the limit of short and strong pulses (e.g. neglecting 
the paths where the photons are emitted within the 
pulse events). Finally, the non-trivial upper and lower 
bounds for the strong square pulses with finite duration 
were obtained. This kind of information is useful in 
experiments, where pulses are neither infinitely short 
not infinitely strong. Our approach can be applied to 
the theoretical study of other types of non-linear spec- 
troscopy such as three level systems, systems undergoing 
stochastic dynamics [3j, or Josephson junction qubits 
[^] controlled by microwave radiation, where the strong 
dependence on the contribution of coherence was already 
experimentally proved. 
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Israel Science Foundation. 
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APPENDIX A 

In this appendix we give the detailed derivation of Eq. Consider the semiclassical trajectories Eq. (pij) of P„ 

from the two pulses . First let us sort the terms according to 

n n— 1 n— 1 

Q— a— a— 

n— 2 ^ n n— 1 

+ E(-l<:*T'^l5)(e|C/£!o)l5) + e-^^ E(.9l<r;i-)(^l<o)l5> + E (-I^I)"' 1^) (^I^^Ko) 1^) 

a=0 {a=0 a=0 

n-1 Ti-2 

- E(5it^ir*r)"^i5)(^it^£lo)i5) ~ E(^i<rr'^i5)(eic^£lo)i5) (60) 

Now we consider the sum of the first two terms 

n n—1 



E(5l<Z^)(5l<o)l5) + E(^l<X)"^l5)(5l<o)lff) 

q:=0 a=0 



+ 



= PnH9\U[Zo)\9) + n!il(.9|f/£o)l3> + • ■ - E ^nia(ff|t^£!o)lff)' (61) 

Similar manipulations with the second two terms lead to : 

n— 1 n— 2 n—1 

E(5|C^£r)"^l5)(^I^So)l.9> + E(^|f^(::rr'^l3>(^|t^£!o)l5) = E ^nii-.(e|<o)l.9) (62) 

a— a— a— 

Combining Eqs. ([CTjl and ((^2) we get : 



E ^.'-.(5l<o)l5) + E ^nii-.(e|<o)lff) 



^nM3l<.o)lff) + {^n-i(e|<,o)lff) + ^n-i(5|C^S,o)lff>} + {n!^-2(e|C/£o)l3) + n!^-2(5|C^£o) 1.9) }+• ' ' = E ^n-aP^ 

a=0 

(63) 

Now we concentrate on the terms multiplied by the factor e ^'^ in Eq. ([50)1 . Let's add and subtract the two following 
terms : 

n n— 1 

E(5l<r)'l.9)(3l<o)l3> + E(-l^fc:r'^l5)(5|t^ffo)l5). 

a=0 a=0 

We get 



e 



^' n n—1 

^ E(5l<:*5l^)(e|<o)l3) + E(^l^£r)"^l^)(-l<o)l5) + E(5l<r?l5)(ffl<o)l5) 



a=0 
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n-l 

a=0 



n-l 



a=0 



a=0 



n— 1 n— 2 



(64) 



a=0 



Q = 



But the first 4 paths of Eq. ([64)) are just the semiclassical part of probabiUty of emitting n photons from the two 
pulses attached together 



pciaj.i. = E(5l<:.S i\9){9\ + |e)(e|)<o)|3) + E(^l<:i""' + 1^) (e|) <o) l^) 



(65) 



a=0 



a=0 



(compare with Eq. (HH])). And the last 4 paths of Eq. are equal to Eq. (|6ip +Eq. (p^ . Putting all this information 
together we obtain 



Q=0 I. k=0 

Finally, by addition and substraction of the coherence trajectories: 



(66) 



a=0 



to pcia,i2ii obtain Eq. ([Ml). Using Eq. 
paths now oscillate in A with loq. 



it should be taken into account however that not all the coherence 
I 



APPENDIX B 

The Rotating Wave Approximation (RWA) [2l| con- 
sists of neglecting the non-resonant processes of rising 
from \g) to |e) by emitting a photon and falling from |e) 
to \g) by absorbing a photon. Switching to the rotating 
frame by applying the transformation ^(t-t',^,;^^) defined 
below, it is possible to suppress any time dependence in 
the Bloch equation Eq. ([T]). As a result the following 
time independent equation is obtained 



and the transformation A(t_ 



^(t) = L + T ai^t) 

where 

t' - is the initial moment. 



L = 



(67) 








-in 

2 


in 

2 








in 


-in 


2 


2 


-in 

2 


in 

2 







in 

\ 2 


-in 

2 





2 ' ' 



(69) 



A 



{t — t',4i,UJL) 



/ 1 
1 

Vo 



is given by 








I 

(70) 

where </> is the phase of the laser at the initial moment t' . 
In the new representation the calculation of the Green 
function is straightforward (see Eq. ([74]) ). Represent- 
ing the solution to the time independent Bloch equation 
Eq. (j67p in the rotating frame as the infinite iterative 
expansion in F we find the following expression for n- 
photon-propagator within RWA 



(68) [/, 



(«-*') 
Hence 



^(t-t„)f •••f^(ii-t')dii---dt„. 

(71) 



(„) 



r(") 



(72) 



where 0{t') 



^(S>,(i>,ujL)'^{t') is the initial condition. Fi- 
nally for obtaining a)^ we have to apply the inverse 
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transformation 



A; 



(") 



A< 



(73) 



From Eq. ([75)1 one easily makes the following conclu- 
sions: (i) the initial state of the molecule is multiplied 



by A(o, 



where 0i is the initial phase of the laser at 



the beginning of the first pulse, which shifts the initial 
coherence phase by —ipi (ii) the delay period propagators 
are now multiplied by ^(o,02,i^l) from the left due to the 
second pulse and by A 



(T,<t>l,u>L) 



from the right due to the 
first pulse {(j)2 is the initial laser phase at the beginning 
of the second pulse). Clearly, this leads only to an 
additional phase shift {Tlol + </>i — 02) of the coherence 
terms. Therefore calculating the photon statistics for 
the square pulses we rewrite the Eqs. (|20)) . (|2T|) . (|22p with 
the following modifications: 

1) In the definition of the n-photon- propagator Eq. ^ 
the Green function defined as the time ordered exponen- 
tial are replaced by 



which are Green functions for the time intervals inside 
the pulses within RWA. 



2) The initial state of the system must be replaced by 

^(O,0,u;i,)Cr(O) 

3) The coherent terms are multiplied by addi- 
tional phase factor e'("^'^+'^i-"^2)^ 

Remark : If in experiments the initial phases of the 
pulses are random variables, it's necessary to replace all 
the phase factors with their ensemble averages. 



Thus summarizing we have 



p„ ^ pCla ^ g-Ar/2 U(A+T)^o+0i-02„Coh ^ ^ ^ 



(75) 



g{t - t') = exp {t - t')L 



(74) where 



n— 1 n— 2 



n— 1 n— 2 

(76) 

and 



a^"'^ = E(5|t^£r) l^)(^|t^(So)l5> + E(^|f^£r)"^l^)(-|t^(£o)l5), (77) 



with \g) = A(o,0,,^^)|.g). 



APPENDIX C 

Here we show the derivation of exact expressions for Pi and P2 for two equal square pulses obtained within RWA . 

For Pi : The one-photon-propagator in (0,t) may be decomposed to the sum of 4 different terms 

1 

7-r(l) -Y^rrK") 7-r(l-«)7-r(0) tjW , rr(0) rr(0) jAa) rr(l-a)] (jc,-. 

^(t,o) - (tMr {HMr (.t^Mr (ti,o) ^ {t^r (HMr {t2,tir [tiA) \ ^'°> 



fe=0 



Inserting the closure relation Eq. (|15p between every two propagators of Eq. ([78)) and applying RWA we obtain the 
probability of emitting a single photon using the trajectories notation 

P?^ = E {(5l?^41l5)(5lCr^l3) + {9\ul:he){e\U^lX\g)e-^^] 



15 



+ {9\U^-l\9){e\ui%\9) (1 - e-^^) + {e\ui^l\e){e\uto%)e-'^^ + {e\U^-l\9) {9^1%) ■ 

«r = j2(9\uil\^){c\uto''^\9) + {e\U^'JtJc){c\ul%). 

After some tedious algebra using Eqs. ([3]), ((69)) . ((7T|) . ([74]) we finally obtain : 

Pi = ai + 6ie"'^ + cie"^ cos [ujq (T + A)] . 

Where 



0.1 



16y 



7(1-^^) 



an + ai2 cosh ( ^ j + ai3 sinh ) + au cosh(Ty) + ais sinh(Tj/) 



Vl -4172, 



(79) 
(80) 

(81) 

(82) 
(83) 



ail =yiy^~ 1) [4 {y^ - 3) + ST {y^ - l)] , ai2 = (r?/^ + 8?/ - T - 16) , ai3 = 4 [-(T + 3)2/"* + (T + A)y^ + 3] , 
ai4 = (T + 4)2/5 + 6T2/3 + (T _ 20)y, ais - 2 [(2r + 5)y^ + 2{T - 5)y^ - 3] . 



^-1 = 1^(1-^^)' 



16y^ 



And 



Cl 



16y 



7(1-^^) 



-3Ty + 2Tycosh ( ^ ) + 12 sinh [ ^ ) + Tj/cosh(T2/) - 6sinh(ry) 



cii + ci2 cosh ( ^ j + ci3 sinh ) + cm cosh(T2/) + cis sinh(Ty) 



(84) 



(85) 



cn - -2?; [2 (y^ - 3) + T (y^ + 2y^ - 3)] , ci2 - 4y [{T + A)y^ - T - 8] , C13 = 4 [Ty^ + (2 - T)y2 _ 6] , 
ci4 = 2y ((T + 2)2/4 _ 82/^ - T + lO) , C15 = 4 [Tz/* - (T + 1)2/^ + 3] . 



For P2 : 

The two-photon-propagator may be decomposed to the sum of 10 terms. 

2 1 



Q=0 /3,ct=0 



( Since u['^\^ = there are only 8 non-zero terms.) This leads to the following expressions for ^2^1^ and ^1^°^ ; 



pf^ - E { 

fc=0 



{9\u\-l^)\9){9\ut.S^\9) + (.9l?^£!..)|e)(e|i/(to?l5>e-"''} + (elulZt^g) {e\u[^l,^^ ( 



1 - e 



-rA\ 



,Coh 



E(.9|t^(St.)l^)(^|t^(to)l5)+E(^l^(S 



(87) 



fc=0 



fc=0 
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Calculating the matrix elements with the help of Mathematica yields: 

P2 = a2 + b2e-'^ + cae"^ cos [ujq [T + A)] 

where 



a2 



64?/ 
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0.21 + 022 cosh ( ^ j + sinh ) +024 cosh(Ty) + 025 sinh(r?/) 



(88) 



(89) 



021 = (9r2 + 40r + 24) / - 2 (gr^ + 56T + 17) y'' + 9 (T^ + BT - 12) + 126, 

022 = (r^ - 32) j/*^ + 4(r + 32)?;'* - (T^ + 20T + 64) - 192, 

023 = 2y [(?/2 _ 1) r2y2 + T (3?;4 - B?;^ - 3) - 4 (B?/^ - 39?;2 + 5l)] , 

024 = (r^ + BT + B) y*^ + 2 (3T2 - 6T - 47) ?/" + (T^ - b2T + 172) + 66, 

025 = y [4 (?/' + 1) T2y2 + (i7j^4 „ 58y2 _ r _ 4 (4^4 ^ 9^2 _ 5^)] ^ 



64?/ 



10 



{y'-^f 



&21 + &22 cosh 



Ty 



623 sinh 



T?/ 



624 cosh(r?/) + &25 sinh(T?/) 



(90) 



621 = 3 (3T2 - B) - IB (r^ - 7) y'' + 9 (T^ - 28) y"^ + 126, 622 = - (T^ - 32) + 2 (T^ - 96) y^ - {T^ - 3B4) y^ 
623 = -6Ty (?/2 _ 1)^ 624 = (T^ - 8) / _ 2 - 33) y^ + (t^ - 132) y^ + 66, 625 = -15Ty (y^ _ i)^ . 



And 



C2 



32y 



10 



y' - 1 



C21 + C22 cosh 



Ty 
2 



C23 sinh 



Ty 



C24 cosh(r?/) + C25 sinh(T?/) 



(91) 



C21 = -r(T + 4)?/^ + 2 (ST^ + 12T - 15) y^ - 3 {ST^ + 12T - 44) y^ 

C22 = - (T^ + 2r - 32) y4 + (T^ + lOT - I6O) y^ + 192, 

C23 = - (T^ - 32) y^ + (T^ + 2r - 156) y^ + 6{T + 34)y, 

C24 = T{T + A)y^ - 2(11T + l)y'* + (-T^ + 26T + 2B) y^ - 66, 

C25 = y [2 (y' - 1) r2y2 + T (-3y4 ~ 4y2 + 15) - 2 (By* - 39y2 + 51) 



126, 
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